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ABSTRACT 



^ , Context. This paper is the last in a series devoted to the analysis of the binary content of the Hipparcos Catalogue. 

Aims. The comparison of the proper motions constructed from positions spanning a short (Hipparcos) or long time (Tycho-2) makes it possible 
t^J- to uncover binaries with periods of the order of or somewhat larger than the short time span (in this case, the 3 yr duration of the Hipparcos 

mission), since the unrecognised orbital motion will then add to the proper motion. 
CN| Methods. A list of candidate proper motion binaries is constructed from a carefully designed \ 2 test evaluating the statistical significance 
^~ H of the difference between the Tycho-2 and Hipparcos proper motions for 103134 stars in common between the two catalogues (excluding 
components of visual systems). Since similar lists of proper-motion binaries have already been constructed, the present paper focuses on the 
evaluation of the detection efficiency of proper-motion binaries, using different kinds of control data (mostly radial velocities). The detection 
, rate for entries from the Ninth Catalogue of Spectroscopic Binary Orbits (S B 9 ) is evaluated, as well as for stars like barium stars, which are 
O |' known to be all binaries, and finally for spectroscopic binaries identified from radial velocity data in the Geneva-Copenhagen survey of F and 
G dwarfs in the solar neighbourhood. 

Results. Proper motion binaries are efficiently detected for systems with parallaxes in excess of ~ 20 mas, and periods in the range 1000 
- 30000 d. The shortest periods in this range (1000 - 2000 d, i.e., once to twice the duration of the Hipparcos mission) may appear only 



o 

, as DMSA/G binaries (accelerated proper motion in the Hipparcos Double and Multiple System Annex). Proper motion binaries detected 
• • ■ among S B 9 systems having periods shorter than about 400 d hint at triple systems, the proper-motion binary involving a component with a 
. ' l° n g er orbital period. A list of 19 candidate triple systems is provided. Binaries suspected of having low-mass (brown-dwarf-like) companions 
, are listed as well. Among the 37 barium stars with parallaxes larger than 5 mas, only 7 exhibit no evidence for duplicity whatsoever (be it 
■ spectroscopic or astrometric). Finally, the fraction of proper-motion binaries shows no significant variation among the various (regular) spectral 
' classes, when due account is taken for the detection biases. 
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1. Introduction 



Among the 118 218 entries in the Hipparcos Catalogue (ESA 
19971) . about 18 000 (collected in the Double and Multiple 
Systems Annex - DMSA) involve binary or multiple stars, 
which break down as follows (Sect. 1.04 of lESAl \l99% : 
13211 entries are component solutions (i.e., visual double 
or multiple systems, DMSA/C), 2622 are acceleration so- 
lutions (i.e., with a variable proper motion, or DMSA/G), 
235 astrometric binaries have a known orbit (DMSA/O), 288 
are Variability-Induced Movers (VIMs, or DMSA/V) and fi- 
nally 1561 are stochastic solutions (DMSA/X), a number of 
which might be binaries. With the exception of component 
solutions, all the other kinds of binary entr ies have been 
re-evaluated in the present series of papers (IPourbaix et al 



Senior Research Associate, F.N.R.S., Belgium 



2003; Platais et al 



Jancart et al. 2005 



2003; Pourbaix 2004; Jorissen et al. 2004 



). Among th e 235 DMSA/O systems, 122 



have been re-investigated by iJancart et al.1 (|2005) in light 



of their spectroscopic orbital elements, and 8 orbital solu- 
tions were rejected b ecause of their low quality (Sect. 4.3 



of Jancart et al. 2005). At the same time, 51 new DMSA/O 



were found using their spectroscopic elements published af- 
ter the completion of the Hipparcos Catalogue, and collected 
in the Ninth Catalog ue of Spectroscopic Binary Orbits (Sgo; 
Pourbaix et all 12004b . For the 188 DMSA/V entries that are 



long-period variables, the re-processing of their Hipparcos data 
(using a variab le chromaticity corr ection) led to the rejection of 
161 of them dPourbaix et al.ll2003l) . 

Even when spectroscopic orbital elements are not avail- 
able, for each star from the Hipparcos Catalogue (with the ex- 
ception of the DMSA/C component solutions), it is possible 
to quantify the likelihood that their Intermediate Astrometric 
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Data (IAD; van Leeuwen & Evanslll998 ) are better fitted with 
an orbital model than with a standard single-star model (i.e., 
a 5 -parameter model, including the parallax w, the a and 
S positions at the 1991.25 reference epoch, and the proper 
motion in a and <5). Hipparcos data are, however, seldom 
precise enough to derive the orbital elements from scratch. 
Therefore, when a spectroscopic orbit is available beforehand, 
it is advantageous to import e,P,T from the spectroscopic 
orbit and to derive the remaining astrometric elements (as 
done by [P ourbaix & Jorissen 2000; Pourbaix & Boffinl 12003; 
Ijancart et al.l 12005 ). If a spectroscopic orbit is not available, 
trial (e,P,T) triplets scanning a regular grid (with 10 < 
P(d) < 5000 imposed by the Hipparcos scanning law and the 
mission duration) may be used. In practice, it is even possible 
to avoid scanning e and T, and to restrict the grid to the or- 
bital period P, assuming circular orbits. This choice greatly re- 
duces the computing time, since test cases revealed that binary 
systems with moderate eccentricities are nevertheles s detected 
even though they are fitted with circular orbits (see [Pourbaixl 
20041 for details). Very eccentric orbits might be missed, of 



course, but their small number does not justify the tremendous 
increase in computer time required to detect them . This metho d 
has been applied to barium stars by lJorissen et al ] (l2004l2005h . 
and has been shown to provide the best detection efficiency if 
w > 5 mas and 100 < P(d) < 4000. 

The present paper is the last one in our series reprocess- 
ing binaries in the Hipparcos Catalogue. It deals with bina- 
ries that are identifiable on the ground of their variable proper 
motion. The DMSA already contains 2622 such DMSA/G 
solutions. Thispaper aims at extending the list of proper- 
motion binaries] by comparing t he Hipparcos pro per motion 

l2000alfbh. follow 



with the Tycho-2 proper motion (H0g et al- — 

ing the idea put forward bvlwielenl (1 19971) and IWielen et al 



( 199 8lll999l).lK arjlan & Ma karovl d2003l) . lMakarovl d2004l) . and 
Maka rov & Kaplan! d2005l) presented further developments and 
applications. The Hipparcos proper motion, being based on 
observations spanning only 3 yrs, may be altered by the or- 
bital motion, especially for systems with periods in the range 
1500 to 30000 d whose orbital motion was not recognised by 
Hipparcos. On the contrary, this effect should average out in 
the Tycho-2 proper motion, which is derived from observa- 
tions covering a much longer time span. This method works 
best when applied to stars with parallaxes in excess of about 
20 mas. The method is fully described in Sect. 

While the p re sent paper was being prepared, 
Makarov & Kapla 2 d2005l) issued the list of proper- 



motion binaries at which we were aiming. Wielen et al 
(1999) published a list of proper-motion binaries (col- 
lected in the DMUBIN database: http://www.ari.uni- 
heidelberg.de/datenbanken/dmubin), based on the comparison 
of Hipparcos proper motions with several other catalogues, as 
well. We nevertheless believe that the present analysis remains 
timely, for the following reasons: (i) the statistical method 
used here to identify proper-motion binaries is different from 



1 The concept of proper-motion binaries, which is the topic of the 
present paper, should not be confused with the more traditional one of 
common proper-motion pairs. 



the one of M akarov & Kaplan ( 2005 ). Our method takes into 
account the correlation between the a and S components of the 
Hipparcos proper motion, as well as between the Hipparcos 
and Tycho-2 proper motions, though on a statistical rather than 
individual manner for the latter, (ii) We illustrate the astrophys- 
ical promises of the method by applying it to specific classes 
of stars like barium stars (Sect. 15.21 ). Am stars, or Wolf-Rayet 
stars (Sect. |6]). Among those classes, binaries are sometimes 
difficult to find using the traditional spectroscopic approach. 
For instance, in the case of barium stars, the periods may be 
quite long, (iii) In Sect. 15.11 we test the d etectability criteria 
of pr oper-motion binaries devised by iMakarov & Kaplan] 
( 2005) by looking for such binaries in the Ninth Catalog ue of 
Spectroscopic Binary Orbits (Sm: iPourbaix et al.ll2004h . (iv) 
The cross-check between the list of candidate proper-motion 
binaries and the corresponding radial-velocity data (available 
for a la rge sample of sta rs from the HIPPARCOS/CORAVEL 
survey: IUdry et al 111 9971) allows us to identify spurious proper- 
motion binaries (due to perspective acceleration) or binaries 
with low-mass (brown-dwarf-like) companions (Sect. l5"T3l >. (v) 
We evaluate the proportion of proper-motion binaries among 
the different spectral classes (Sect. [6]) and do not confirm 
Makarov & Kaplan's (2005) suggestion that the fraction of 
long-period binaries is larger among solar-type stars than 
among the other spectral classes. 

2. Statistical method 



The Tycho-2 Catalogue dH<3g etal.ll2000allbT) is an extension of 
the Tycho Catalogu e collected by the sky mapper of the ESA 
Hipparcos satellite dESAlll997l) . It contains astrometric posi- 
tions and proper motions for 2.5 10 6 stars. The proper motions 
were extracted from the Hipparcos and Tycho-2 catalogues 
for those stars with only one component, thus excluding the 
DMSA/C and DMSA/O entries, as well as Hipparcos objects 
that have double entries in the Tycho-2 catalogue. We ended up 
with 103 134 stars. The difference between the Hipparcos and 
Tycho-2 proper motions is then evaluated through a \ 2 test: 



X' 



= x l s- 1 



X. 



(1) 



where X is the two-dimensional vector containing the differ- 
ence between the Hipparcos and Tycho-2 proper motions, the 
first component corresponding to the right ascension a and the 
second to the declination 5. The 2x2 covariance matrix S on 
the proper-motion differences is obtained from the 4x4 covari- 
ance matrix £ on the proper motions themselves by the relation 



S = A £ A\ 



where A is a 2 x 4 matrix such that 



AY 



10-1 
1 0-1 



(2) 



(3) 



where Y is the 4-dimensional vector containing the a and S 
components of the Hipparcos proper motion followed by the 
a and 8 components of the Tycho-2 proper motion. The co- 
variance matrix £ associated with the proper-motion vector Y 



A. Frankowski et al.: Proper-motion binaries in the Hipparcos Catalogue 



3 



/ a aH PH a aH &8H PHT,u &aH & aT \ 

PH &aH OSH &SH PHT.8 OSH V ST 

PHT.a &aH &aT f7^ T px <J aT &5T 

\ PHT.S C5H &ST PT CT q t &8T &ST / 



is expressed by Eq. |U Here a a n and asn are the standard er- 
rors of the Hipparcos proper-motion a and S components, from 
fields H17 and HI 8 of the Hipparcos Catalogue, and a a T and 
<TgT are the respective standard errors of the Tycho-2 proper 
motion, from fields T17 and T18 of the Tycho-2 Catalogue. 
The correlation coefficient p# between the a and S compo- 
nent of the Hipparcos proper motion is taken from field H28 of 
the Hipparcos Catalogue. The analogous correlation pt for the 
Tycho-2 catalogue is supposed to be null, given the different 
approach used to construct the Tycho-2 catalogue. The correla- 
tion coefficients pht,o and pht,s between the Hippa rcos an d 
Tycho-2 proper motions were evaluated by Hei g et alJ (I2000ah . 
based on a sample of supposedly single stars. To quantify these 
correlations, these authors define a parameter R (separately for 
the a and 6 components; when no ambiguity is possible, the 
subscript will therefore be omitted, for the sake of simplicity) 
as 



R 2 = 



'H 



and 



~2 2 

°d = a H 



<J T — 2pHT &H &T + f c 



(5) 



(6) 



a c being a 'cosmic error' . The quantity a c in fact accounts for 
the presence of unrecognised astrometric binaries, which arti- 
ficially increases <Jd when comparing proper motions from two 
different catalogues dWielenlll997llwielen et al.lll998h . This is 
precisely the effect we are after. Therefore, including the cos- 
mic error in the evaluation of R, which is itself used in the 
X 2 test (Eq. Q]), would in fact introduce a bias in the method. 
Therefore, a new quantity r has been evaluated, closely follow- 
ing the definition of R (Eq. [5]), but without the cosmic error 
term. It is listed in Table Q] an d is de rived from the R values 
listed in Table 1 of lH0g et a D (EoOOa) using the following ex- 
pression: 

^2 ! 

-■ (7) 



H 



1 

w 



The parameter R was obtained by lH0g et al. | d2000al) for 
bins of 1 magnitude from Vt — 7.0 to 12.0 (restricted to sub- 
samples of non-binary stars), by computing ad and the rms 
values of the internal errors an and ctt of the proper motions 
given in the two catalogues, plus an estimate of the cosmic er- 
ror. It is important to stress at this point that no R factor is 
available for the brightest stars (Vt < 7). These bright stars 
were not considered in the derivation of R as the vast majority 
of them, having accurate Hipparcos data, determine the sys- 
tems of the other catalogues, from which the Tycho-2 proper 
motions are derived. 

The parameter r thus expresses the ratio between the stan- 
dard deviation ad of the difference in Hipparcos and Tycho-2 
proper motions and the root-mean square of the internal errors 
for these two catalogues. In the absence of any correlation, r 



Table 1. The quantities r a and rs as a function of the Tycho 
magnitude Vt (see Eq. [7} 



V T 


r a 


rs 


< 7.0 


1. 


1. 


7.0-8 


0.90 


0.95 


8.0-9 


0.92 


0.94 


9.0- 10 


0.89 


0.93 


10.0- 11 


0.95 


0.99 


11.0- 12 


1.03 


1.09 


> 12.0 


0.98 


1.13 



should be unity; in the presence of correlation, p^T ma y be 
expressed 



PHT &H or 



2 \ I J 1 



4). 



(8) 



It should be noted, however, that a value of r different from 
unity may also indicate that the 'internal errors' an or cry 
for the proper motions from the Hipparcos and Tycho-2 cat- 
alogues, respectively, were incorrectly evaluated. Equation [7] 
may indeed be rewritten a\ ~ r 2 a^ + r 2 a\. The two inter- 
pretations of r are in fact algebraically equivalent, as it will 
become apparent in Eq.|9] It is shown in Fig. Q] that the inclu- 
sion of r is necessary to ensure that the \ 2 variable of Eq. [Tire- 
ally follows a x 2 distribution. The correctness of our statistical 
treatment is also illustrated in Sect. [6] in relation with Fig.fTTI 
where the rate of false detections is exactly the one expected. 

After some algebra, it is then possible to show that 



r a ( ff „H + °qt) Ph r a r s a aH o&H 
PH r a r s a aH a SH rj {a 2 SH + a 2 ST ) 



(9) 



It must be noted that in the derivation of Eq. [9] to ensure that 
the determinant of the S matrix is always positive (ensuring the 
positiveness of the x 2 as defined in Eq.Q}, the correlation term 
Ph &aH &8H in Eq.|]must be replaced by p H r a rs a aH asn- 
This modification follows naturally from the realisation that the 
r a and rg factors account for the possible mis-evaluation of the 
internal errors a a n, 

3. Results: Candidate proper-motion binaries 

The confidence level of 0.9999 on the x 2 probability to flag a 
star as binary has been set to minimise the fraction of false de- 
tections. This threshold translates into x 2 = 18.4206 (with two 
degrees of freedom). The insert in Fig.[TJreveals that 1734 stars 
out of 103 134 have Prob(x 2 ) > 0.999^1 representing 1.68% 
of the total sample, whereas the expected number of false de- 
tections is 0.01% of 103 134, or only 10 stars. In the remainder 



2 In the remainder of this paper, Prob(x 2 ) denotes the left-sided 
cumulative % 2 probability, i.e., from to x 2 - 
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0.4 0.6 

Prob(x 2 ) 

Fig. 1. Comparison of the cumulative frequency distribution of 
the x 2 values from Eq. [Devaluated with (thick solid line) and 
without (thin solid line) the correction term r from Eq.|7] Note 
how well the ? -corrected distribution follows the diagonal (dot- 
ted line), as it should in the absence of binaries. The brightest 
stars with Vr < 7 were not corrected, and correspond to the 
dashed line which deviates markedly from the diagonal. The 
insert is a zoom on the upper right portion of the diagram, for 
the complete sample of 103 134 stars. 

of this paper, stars fulfilling the condition Prob(x 2 ) > 0.9999 
will be called (candidate) 'A/i binaries'. Among these 1734 
candidate binaries, 729 were already flagged in the Hipparcos 
DMSA as exhibiting a non-uniform proper motion (i.e., 162 
DMSA/G solutions with 9 parameters: the 5 usual astrometric 
parameters, plus 4 parameters describing the first and second 
derivatives of the proper motion, and 567 DMSA/G solutions 
with 7 parameters, corresponding to a linear variation of the 
proper motion). 

Several more proper-motion binaries actually exist in the 
Hipparcos Catalogue, since 1831 among the 2560 DMSA/G 
stars present in our sample were not flagged as binaries by the 
X 2 test. This is in part a consequence of the fact that the proper 
motion provided by the Hipparcos Catalogue for DMSA/G 
stars already includes some correction for the orbital motion, 
so that the difference between the long-term (i.e. Tycho-2) and 
short-term (i.e. Hipparcos) proper motions is not necessarily 
significant. In total, we thus identify 3565 candidate proper- 
motion binaries (=1831 DMSA/G stars not reaching the 0.9999 
confidence level + 1734 stars above that level; among the lat- 
ter, there are 1734 - 729 = 1005 candidate proper- motion bina- 
ries not already flagged as DMSA/G). These 3565 objects are 
listed in Table [2] only available in electronic form at the CDS, 
Strasbourg. This number is close - albeit not identical - to the 
3783 candidat e binaries (among which 1 161 are not DMSA/G) 
identified bv lMakarov & Kaplan! (12005b in their Tables 1 and 2. 



Before comparing our results with those of Ma karov & Kaplan] 
(2005), it should be mentioned that among DMSA/X (stochas- 
tic) solutions, the fraction of candidate proper-motion binaries 
(45/1403 or 3.2%) is only slightly larger than among the full 
sample (1.68%), suggesting that long-period binaries are not 
the major cause of stochasticity. 

4. Comparison with Makarov & Kaplan results 

Mak arov & Kaplan! d2005 ) adopted a very simple criterion to 
flag a star as a candidate proper-motion binary, namely that 
the difference between the Hipparcos and Tycho-2 proper mo- 
tions exceeds 3.5 times their root-mean square error for at 
least one of the two components. For a Gaussian distribution, 
this threshold corresponds to a fraction of false detections of 
2 x 4.7 10~ 4 - 10~ 3 , or 100 stars in a sample of 10 5 . Our 
threshold is a factor of 10 more stringent. This is immedi- 
ately apparent in Fig. [2] which compares the x 2 values for 
Makarov & Kaplan's candidate proper-motion binaries (shaded 
histogram) with our 0.9999 confidence level corresponding to 
X 2 > 18.4206. A 0.999 confidence level (corresponding to 
the false detection rate of Makarov & Kaplan) translates into 
X 2 > 13.8155 (for two degrees of freedom), which indeed 
matches the left edge of Makarov & Kaplan's x 2 distribution 
of their candidate proper motion binaries. 

It is in the range 12 < x 2 < 25 that the two sets dif- 
fer, as may in fact be expected (Fig. |2): in that range, 316 
stars were flagged as proper-motion binaries by Makarov & 
Kaplan, but not by us, while 138 were flagged by us and not 
by Makarov & Kaplan. This difference mainly results from 
the different criteria used in these two studies: if the differ- 
ence between the Hipparcos and Tycho proper motions is about 
3.5c for one component only, the corresponding x 2 will on ly 
amount to 3.5 2 = 12.25, well below the detection thresh- 
old used in the present work. Conversely, the difference in 
each component might be slightly smaller than 3.5cr, but the 
combined x 2 (Eqs. Q] and |9]l may be larger than 18.4 (for 
instance 3.1 2 + 3.1 2 = 19.2). A good example thereof is 
HIP 1 15929 not flagged by Makarov & Kaplan, despite its x 2 
of 26.12 yielding Prob(x 2 ) = 0.9999978 in our method. With 
/i a *(HIP) = -5.29 ± 1.08 mas y" 1 and /i Q * (Tycho — 2) = 
— 10.2 ± 1.2 mas y^ 1 , the difference in fj, a * amounts to 
|A^|/[e 2 (HIP) + e 2 (Tycho- 2)] 1 / 2 = 3.0, and with 
Ata(HIP) = 36.52 ± 0.81 mas y~\ ^ (Tycho - 2) = 32.1 ± 
1.1 mas y _1 , the difference in 6 amounts to | A/15 |/[e 2 (HIP) + 



e 2 (Tycho - 2)] 



1/2 



3.2. An additional factor acting in the 



same direction is that, with a Tycho-2 magnitude Vr — 9.01, 
HIP 1 15929 belongs to the magnitude range where the r a and 
rs parameters (Table Q]) have the largest impact on the x 2 
value, increasing it by about 20%. Incidentally, HIP 115929 
is also flagged as a A/i binary in the DMUBIN database 
JWielenetalJll999l) . 



5. Comparison with other samples of known 
binaries 

Our list of candidate proper-motion binaries identified in 
Sect. [3] can be confronted with several independent data 
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Table 2. The first 20 entries in the list of 3565 DMSA/G or A/i binaries at the 0.9999 confidence level. The full list is only 
available electronically at the CDS, Strasbourg. \ 2 an d Prob(x 2 ) characterise the significance of the difference between the 
Hipparcos and Tycho-2 proper motions (see Eq. [TJ. The 'DMSA' column refers to the Double and Multiple Systems Annex 
of the Hipparcos Catalogue (ESA 1997): '5' stands for single-star (5-parameter) solution, '7' and '9' for accelerated proper 
motions (7- and 9-parameter solutions), and 'X' stands for stochastic solution. The columns labelled 'S pectrosc. bin.' provid e 
information about the possible spectroscopic -binary nature of the star, using da ta from the Sr9 catalog ue (iPourbaix et alJl2004h . 
th e Geneva-Copenhage n survey of F and G dwarfs in the solar neighbourhood (INordstrom et al.l l2004). or the KIII-MIII sample 
of lFamaey etaD |2005). In those columns, 'B' or 'b' means that the star is flagged as a spectroscopic binary in these catalogues 
(with Prob[x 2 (K)] < 0.001 and 0.001 < Prob[x 2 (X)] < 0.01, respectively), 'O' that a spectroscopic orbit is available, 'c' 
that it is present in the catalogue but that it is not flagged as binary. The column 'Remark' provides a comment relating to 
specific properties of the star: long-period variable (LPV); barium star (Ba); triple system (triple) or candidate triple system 
(triple?); suspected binary with low-mass companion (sB-lm?); spurious proper-motion binary due to perspective acceleration 
(perspective ace). 



HIP x 2 ProRx 3 ) DMSA Spectrosc. bin. Remark 
S B 9 Gen-Cop K-M " 



62 


3.70 


0.84276 


7 








68 


573.65 


1.00000 


7 








93 


32.19 


1.00000 


5 




B 




104 


9.48 


0.99126 


7 








171 


1695.38 


1.00000 


X 


O 


B 




290 


91.51 


1.00000 


5 




B 




305 


42.75 


1.00000 


7 




c 




329 


37.94 


1.00000 


5 








340 


8.28 


0.98408 


9 








356 


43.03 


1.00000 


5 




B 




359 


105.95 


1.00000 


5 








493 


3.74 


0.84589 


7 




c 




611 


11.68 


0.99708 


9 








626 


21.38 


0.99998 


5 








646 


2.98 


0.77418 


7 








648 


5.74 


0.94337 


9 








695 


5.10 


0.92177 


7 








727 


43.81 


1.00000 


5 








741 


13.24 


0.99867 


7 








765 


19.10 


0.99993 


X 









sets available, namely (i) the spectroscopic binaries from 
the Ninth Catalogue of Spectroscopic Binary Orb its (5^9 ; 



Pourba ix et alj 120041) . (ii) barium stars listed by iLii et al 
dl983l) . (iii) the radial-velocity standard deviation available 
from the Geneva-Copenhagen survey of F and G d warfs 



in the solar neighbourhood ( INordstrom et alj 120041) . and 



(iv) the radia l -veloc ity standard deviation available from the 
Famaev et al.1 d2005l) catalogue of K and M giants in the so- 
lar neighbourhood. The latter two data sets are subsamples 
resulting from the major effort to collect radial ve locities for 
about 45 000 Hipparcos 'survey' stars later than F5 dUdrv et al 



1997), using the CORAVEL spectrovelocimeter dBaranne et al 



19791) . These data will, among other things, enable us to check 
whether the proper-motion bina ries are indeed found in t he 
orbital-period range predicted by Ma karov & Kaplan] ( 2005 ) - 
an important observational test for the method, not attempted 
before. As a general remark applying to the remainder of this 
paper, we stress that the orbital inclination somewhat blurs the 
comparison between proper-motion and spectroscopic diagnos- 
tics, since proper-motion binaries seen close to pole-on will not 
be detected as spectroscopic binaries. 



5. 1. S b q entries 



At the time of our q uery, the Ninth Catalo gue of Spectroscopic 
Binary Orbits (Sb9 ; IPourbaix et al.l2004h holds orbits for 2386 
systems, among which 1320 have a Hipparcos and Tycho-2 en- 
try. The period distribution for those stars is presented in Fig.[3j 
along with the frequency of stars for which the Hipparcos and 
Tycho-2 proper motions differ at confidence levels of 0.9999 
and 0.99. Careful inspection of the data shows that the lower 
limit on period for efficient A/i binary detection is located at 
1500 d (see also Fig.|4]i. Proper-motion binaries with periods up 
to 30000 d have been detected, corresponding to the longest 
periods present in the Sgs catalogue. Among the long-period 
proper-motion binaries detected, all eccentricities are equally 
well represented. Figure [4] presents an eccentricity -period di- 
agram for 5g9 stars with parallaxes larger than 10 mas. Note 
that almost all such objects with P > 1000 d have been flagged 
as either DMSA/G or A/i binaries. On the other hand, A/i or 
DMSA/G binaries with orbital periods shorter than about 400 d 
are (candidate) triple systems (see the discussion below and 
Tables[3]and0. 
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a= 0.99 

a= 0.9999 



100 1000 
P(d) 



Fig. 2. Tail of the x 2 distribution containing the proper- 
motion binaries, located to the right of the vertical dashed 
line, corresponding to the threshold used in the present work 
[Prob(x 2 ) > 0.9999, translating into \ 2 > 18-4 for two de- 
grees of freedom] . The can didate proper-motion binaries from 
Makarov & Kaplan! (120051) are identified by the shaded his- 
togram. They are mostly located to the right of the dotted line, 
corresponding to Prob(x 2 ) > 0.999, or \ 2 > 13-8 (for two 
degrees of freedom). The solid line is the theoretical \ 2 distri- 
bution with two degrees of freedom expected in the case of a 
pure single-star sample. 



Makarov & Kaplan (2005) have presented a detailed anal- 



ysis of the relationship between the system's orbital period and 
eccentricity, and the probability that such a system be either 
detected as a proper-motion binary (i.e., with A/i significantly 
different from 0), or as DMSA/G (i.e., with an acceleration 
term in the proper motion, /i). They argued that systems for 
which the latter effect is stronger will be mainly found in the 
period range 3 to 6 yrs (1000 to 2000 d), whereas systems de- 
tectable as A/i binaries, but not as /i binaries will be mainly 
found among binaries with even longer periods. This is exactly 
what is seen in Fig. [4] Systems solely detected as DMSA/G 
and not as proper-motion binaries are mainly restricted to the 
short-period side of the P > 1 000 d range. 



As shown theoretically bv lMakarov & K aplan (2005), sys- 
tems ought to have a parallax large enough to be detected as 
proper-motion binaries. Figure|5]illustrates the incidence of the 
parallax on the detection efficiency: for parallaxes between 10 
and 25 mas, about 50% of the S^s stars with periods larger 
than 1500 d are detected, and that fraction becomes larger than 
80% above 25 mas. Among the 69 systems with P > 1500 d 
and parallaxes in excess of 10 mas, 28 remain undetected, ei- 
ther as DMSA/G or as proper-motion binaries. A closer look at 
those systems reveals however that 15 of them have confidence 
levels in excess of 0.9, and may in fact be considered as proper- 



Fig. 3. Upper panel: Orbital-period distribution for stars from 
the 5g9 catalogue. Lower panel: Fraction of stars in a given 
bin for which the Hipparcos and Tycho-2 proper motions differ 
at confidence levels of 0.9999 (dotted line) and 0.99 (dashed 
line). 



0.8 



0.6 



0,'" 



0.2 



,Prob(x a ) > 0.9999 & DMSA/G 
_ QProb(x 2 ) g 0.9999 & DMSA/G 

,ProbCrt > 0.9»9» & non DMSA/G- 
0.99 < Prob(x 2 ) ^ 0.9999 



0. • ■ ■ 
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Fig. 4. Eccentricity -period diagram for stars from the S^g cat- 
alogue with parallaxes larger than 10 mas. Symbols are as fol- 
lows: large filled symbols: proper-motion binaries detected at 
the 0.9999 confidence level (and flagged as DMSA/G: large 
filled circles; not flagged as DMSA/G: large filled squares); 
open circles: systems flagged as DMSA/G but not fulfilling 
the 0.9999 confidence level; crosses: systems not flagged as 
DMSA/G, lying in the range 0.99 < Prob(x 2 ) < 0.9999. 
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10 20 
parallax(mas) 

Fig. 5. Upper panel: Parallax distribution of Sgo stars with 
orbital periods larger than 1500 d. Lower panel: Fraction of 
those Sg9 stars detected as proper-motion binaries at confi- 
dence levels a of 0.99 and 0.9999. For parallaxes between 10 
and 25 mas, about 50% of the Sgg stars with periods larger than 
1500 d are detected, and that fraction becomes larger than 80% 
above 25 mas. 



motion binaries. There is no reason to expect that all these 69 
binaries should be detected by the A/i method, as at 10 mas 
the detection efficiency of the method is still below maximum 
(as shown in Fig. [5] and further discussed in Sect. |6]on a much 
larger sample). Nevertheless, we note that among the 13 re- 
maining undetected cases, 8 might be explained by one of the 
following reasons: 

- The spectroscopic orbit is of poor quality (HIP 108473S; 

- The system has two visible spectra (SB2 or composite 
spectrum), so that the photocentre is not very different 
from the centre of mass and no orbital motion can be de- 
tected (HIP 2563, HIP 23016, HIP 73449, HIP 114576, 
HIP 111528); 

- The system is a triple or multiple (HIP 20430, HIP 1 14576); 

- The system is eccentric and was observed by Hipparcos 
at a slowly varying orb ital phase (HIP 106786; 
iGontcharov & Kivaevall2002l) . 

As opposed to the situation just considered (long-period 
binaries not detected as A/i binaries), systems with P < 
400 d should not have been detected as proper-motion bina- 
ries. Unless, that is, the short-period spectroscopic orbit is in 



It must be mentioned that HIP 108473, not flagged as a A/i binary 
when comparing Hipparcos and Tycho-2 proper motions, is neverthe- 
less flagged as such when comparing Hipparcos prop er motions to the 
ones computed fro m the Boss General Catalogue l lBoss et ail 1 19371 : 
IWielenetal.lll999l) . 



fact the inner orbit of a wider pair involving a third compo- 
nent. A systematic search for a possible third star in these 
systems has been performed in the Multiple Star Catalogue 
( Tokovinfnll 1997 ) and in the F ourth Catalog of Inte rferometric 
Measurements of Binary Stars (IHarfkopf et al .1200 ll) . which in- 
cludes all published measurements of binary and multiple star 
systems obtained by high-resolution techniques (speckle inter- 
ferometry, photoelectric occultation timing, etc.), as well as 
negative examinations for duplicity. Table [3] lists the 16 short- 
period (P < 400 d) DMSA/G or A/i systems with firm ev- 
idence for being triple or quadruple hierarchical systems. It 
is remarkable that about half of the systems (more precisely 
16/39 or 41%) that are suspected of hosting a third star from 
the proper-motion analysis indeed turn out to be triple after 
close scrutiny. The case of HIP 35487 (R CMa) is exemplary in 
that respect: a detailed analysis of a century-long eclipse tim- 
ing of this system indeed revealed the presence of a third star 
jRibas et alJlioO^ . 

Another remarkable case is HIP 72939, flagged as a 
DMSA/G system but with a period of 3. 55 d listed in ,S B 9 . 
After the present analysis was completed, iFekel et al.1 §006) 
obtained a spectroscopic orbit for a third, outer companion, 
with a period of 1641 d. The revised Hipparcos proper motion 
is now within 1 a of the Tycho-2 one. 

This proportion of 41% triple systems among short-period 
5gs binaries flagged as DMSA/G or A/i binar ies is much larger 



than e xpected from random selection, since Hal bwachs et al 



(2003) find on average 1 triple system for 25-30 systems (3 
to 4%) with F7-K primaries and P < 10 yr. It therefore 
seems reasonable to suggest that the remaining 23 short-period 
(P < 400 d) DMSA/G or A/i systems will turn out to be triple 
systems as well, despite the current lack of evidence. They are 
listed in Table H] A few stars should perhaps be removed from 
this list of candidate triple systems: HIP 21673, HIP 52444, 
HIP 68523, and HIP 7791 1, whose orbital elements are in fact 
considered doubtful (hence, they might not be short-period sys- 
tems after all). 

HIP 77911 deserves further discussion, as it is an outlier 
in the eccentricity - period diagram , with P = 1.26 d and 



0.61 ± 0.2 (ILevato et al.l ll987. who consider their orbit 



as very marginal, however) . A visible companion was found by 
Kouwenh oven et al. I d20Q5h using the ADONIS adaptive optics 



system on the ESO 3.6 m telescope, but it is about 8" away, and 
is thus unlikely to be responsible for the DMSA/G nature of the 
primary star. Either the orbital period of HIP 779 1 1 is in error, 
or the outlying position of HIP 7791 1 in the eccentricity - pe- 
riod diagram r esults from eccentricity pumping by a third star 
in the system dMazeh & ShaharrJl979tlMazehlll990l) . the same 
which is in fact responsible for the accelerated proper motion. 

Considering the 16 firm and the 19 suspected triple sys- 
tems, the frequency of triple systems amount to 35/1320 (or 
2.6%) among the S^s systems considered in the present anal- 
ysis. This frequenc y is consistent with the f requency of triple 
systems reported by lHalbwachs et al.1 (120031) . To conclude this 
discussion, it is interesting to remark that several of the can- 
didate triple stars listed in Table |4] are so-called "twins", 
i.e., short-period (P < 10 d) SB2 systems with components 
of (almost) identical spectral types. It has been argued that 
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many twins are in fact members of higher-multiplicity systems 
(Toko vininl2000i see however Halbwachs et al. (2003) who do 
not share this opinion). 

5.2. Barium stars 



Barium stars represent a very interesting control sample for 
the binary-detection method discussed in this paper, since they 
are kn own to be al l binaries with periods in the range 80 - 
10 4 d (|jorissenll2003l) . Therefore, the detection rate of A/i bi- 
naries among barium sta rs provides a h int on the detection 
biases. The catalogue of iLii et alj d 1983b contains 163 bona 
fide barium stars with an HIP entry (excluding the supergiants 
HD 65699 and HD 206778 = e Peg). This number drops to 
156 after removing DMSA/C and DMSA/O entries. Only 9 
A/i binaries are detected among barium stars at the confi- 
dence level 0.9999 (HIP 13055, 58948, 65535, 71058, 97613, 
98431, 104732, 105294, 112821), 5 of which being moreover 
flagged as DMSA/G. This rather small number of detections 
may appear surprising at first sight, since astrophy sical consid 



erations require all barium stars to be bi naries ( McClure et al 
1980l lMcClurel[l98l iMcClure & Woodsworthlll990l: Ijorissenl 



2003). All but 12 barium stars have parallaxes smaller than 
10 mas, though, which makes the detection of A/i binaries 
among barium stars very inefficient, as discussed in Sects. 15.11 
and 15.31 Among these 12 barium stars with parallaxes larger 
than 10 mas, 6 are indeed detected as A/i binaries at the 0.9999 
confidence level. 

The present astrometric method complements several pre- 
vious studies devoted to the detection of binaries among bar- 
ium stars, using various astrometric and spectroscopic meth- 
ods. Since these studies are scattered in the literature, it is use- 
ful to collect their results in one single table (Table|5]) and figure 
(Fig. |6]). Because the efficiency of the astrometric methods to 
detect binaries depends on the parallax, Table [5] has been split 
in two parts, for barium stars with parallaxes larger or smaller 
than 5 mas. This 5 mas threshold contrasts with the 20 mas 
threshold adopted in Sect. [6] A lower threshold for barium stars 
is motivated by their scarcity at larger parallaxes (Fig. [7]). To 
reach a detection efficiency of proper-motion binaries of the 
order of 30 to 50% around 5 mas we relax the confidence level 
to 0.9 (see the discussion relative to Fig. Q~T|in Sect. |6l com- 
pare with Fig. |7). Objects flagged as A/i binaries by this less 
stringent criterion have a separate flag in Table [5] 

The flags listed in Table |3J referring to various astrometric 
methods, deserve more detailed explanations: 

- Flag 2: The presence of an orbital motion in the astro- 
metric signal may be checked by using trial orbital peri- 
ods a nd assuming a circular orbit (for details, see Sect. [T] 
and |Pourbaixll2004l) . even though the spectroscopic orbit 
is not availabl e. This method has been applied to barium 
stars by iJorissen et al ] (l2004l2005h . and works best when 
m > 5 mas and 100 < P(d) < 4000 (see open triangles in 
Fig.©; 

- Flag 3: An astrometric orbit constrained by the 
spectroscopic orbital elements could be derived 
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Fig. 6. Barium stars with known orbital periods and parallaxes. 
Objects flagged as A/i binaries are represented by filled squares 
ifProb(x 2 ) > 0.9999 (flag 4 in Tabled), and by large filled cir- 
cles if 0.9 < Prob(x 2 ) < 0.9999 (flag 5). DMSA/G binaries 
are represented by large open circles. Large open triangles cor- 
respond to binaries where an astrometric orbit could be derived 
(flag 3) or at least an astrometric signature of binarity has been 
detected (flag 2). Arrows denote stars with only a lower limit 
available on the period. Small filled circles are barium stars that 
were not flagged as binaries by any of the astrometric methods 
to detect duplicity. 



( Pourbaix & Jorissenl 20001 ; Jancart et al. 2005 : open 
triangles in Fig . [6j» ; 

- Flags 4, 5: A/i binaries, most efficiently detected in the pe- 
riod range 1500 - 10000 d. Flag 4 - detection with a confi- 
dence level of 0.9999 (filled squares in Fig.©; flag 5 - with 
a confidence of 0.9 (large filled circles in Fig.©; 

- Flag 6: Binaries detected from their accelerated proper mo- 
tion (flagged as DMSA/G in the Hipparcos Catalogue) cor- 
respond to periods in the range of one to a few 10 3 d (large 
open circles in Fig. [6). 

For the sake o f completeness, spectroscopic binaries (from 
Jorissen et al 1998) are mentioned as well in Table[5]with flag 
1 (flag in the case that no radial velocity variations have been 
detected by CORAVEL monitoring). No flag means that the 
star was not detected as binary by any of the astrometric meth- 
ods and it was not included in the CORAVEL survey. Figure [6] 
is remarkable in demonstrating quite clearly the ability of the 
various astrometric methods to detect the orbital motion all the 
way from 500 d to more than 5000 d, provided that the parallax 
be larger than about 5 mas. The present analysis thus completes 
the series of papers djorissen et al.lll998t IPourbaix & Jorissenl 
20001 IJorissen et al. 1120041120051) devoted to the detection of bi- 
naries among barium stars. Among the 37 barium stars with 
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Table 3. S B a systems with P < 400 d, flagged as DMS A/G or Afi systems, and with firm evidence for being triple or quadruple. 



HIP DMSA/G 



P(d) 



Remark 

SB2: triple system [ Goec king et aljfT 994" Kim et al. 2003) 
SB2 resolved; suspected triple as the systemic velocity seems to 
be variable 

SB2 (F1V + G8IV-V); tripl e system detected by light-travel 
time effect dRibas et alj2002h 



24156 
30630 

35487 



42542 
49136 



60129 
72939 
84886 
86187 



0.42 
6.99 



0.0 
0.15 



1.13 0.0 



3.4 10 4 
48.72 
0.28 

0.81 
71.9 
4792 
3.55 

21.82 
2909 
2.50 



0.49 
0.11 
0.0 

0.0 
0.34 
0.079 
0.0 

0.1126 

0.06 

0.06 



Triple system faarfkopf et al.ll200ll : |Pourbaix et alj[ 2004) 
Quadruple system: two close pairs separated by 0. 147 arcsec on 

the sky, 

in an orbit of period 7 10 ddTokoyininll 19971) 



Triple system (Hartko pf et al 



19921) 



Triple system dMavor & Mazehl 1 19871 : iTokovininl 1 1991 

iFekel et alj20o|) 

Triple system jGriffin & Comelllll997l : lTokovinirj|l997l) 

Triple system: third star found at abo ut 0.3 arcsec; systemic 
velocity of the inne r pair is variable ( Carquillat e t al.l 1 19761 : 
iHartkopf etalfeoOll) 



87655 


7 


0.80 





Triple system: the third star is in an orbit of period 8.9 yrs 
around the close pair (Soderhjelm 1999l:lTokovininlll997h 


88848 


X 


1.81 


0.0 


Triple svstem dFekel et alj|2005t) 




X 


2092.2 


0.765 




90312 


5 


8.80 





SB2 (G2IV + K2IV); suspected non single in Hipparcos 


98954 


9 


198.7 


0.05 


A third component at 2.3" ITokovinin 1997) 


103987 


9 


377.82 


0.086 


Triple system (third star at 0.132" with an estimated period of 
14 vrs; llVIason et alj|200lb 


105091 


7 


225.44 


0.23 


Triple system 




7 


5.41 


0.11 




117712 


5 


7.75 


0.0 


Twin SB2: K3V + K3V; triple system: third star at 4.6" in an 



parallaxes larger than 5 mas, there are only 7 with no evidence 
of orbital motion. 



5.3. Geneva-Copenhagen survey of F and G dwarfs 
in the solar neighbourhood 



The Geneva-Copenhagen catalogue (Nordstro m et al.l |2004) 



contains about 16 700 F and G dwarfs of the solar neighbour- 
hood. For 9431 of those, members of visual pairs excluded, 
CORAVEL radial velocities as well as Hipparcos and Tycho-2 
proper motions are available. It is then possible to confront the 
duplicity diagnostic based on proper-motion with that based 
on radial velocities. The latter criterion is based on the com- 
parison between an internal error estimate e and an external 
error estimate provided by the standard deviation a(V r ) of 
the radial-velocity measurements. From these and the num- 
ber of observations N(V r ), the probability Prob(xy ) that the 
observed radial-velocity scatter is due to measurement errors 
alone may be computed from -y\, = [N(V r ) — 1] ((T{V r )/e) 2 
djorissen&Mayorill988l) . We adopt 1 - Prob(xy r ) > 0.999 
as a threshold for flagging a star with variable radial velocity, 
most likely due to duplicity. Note that, with the above defini- 
tions, Prob(xy r ) is the probability for an object of being single, 



while the proper-motion related Prob(x 2 ) is the probability for 
an object of being a binary. 

As revealed by Fig. [8] the largest fraction (about 15%, rais- 
ing to 22% when adding DMSA/G systems) of proper-motion 
binaries indeed arises among the stars flagged as radial-velocity 
binaries at the 99.9% confidence level. These fractions become 
respectively 29.9 ± 2.4% (= 109/364) and 35.7 ± 2.5% (= 
130/364) when eliminating those stars with parallaxes smaller 
than 20 mas among which proper-motion binaries are not ef- 
ficiently detected. This fraction of only 35.7% proper-motion 
binaries among spectroscopic binaries reflects the fact that 
spectroscopic and proper-motion binaries span different pe- 
riod ranges: Ppm > 1000 d for proper-motion binaries in 
general, and Pa^ > 1500 d for A/i binaries (see Figs. [3] 
and |U). Let Psb be the maximum period of spectroscopic 
binaries detectable with the radial-velocity coverage of the 
Geneva-Copenhagen catalogue, to be fixed later on. The rel- 
ative fractions of those binaries may be computed from the pe- 
riod distributi on of binaries, which has been shown to be log 
- normal fouauermov & Mavorlll99ll : lHalbwachs et al. 2003 
Eggenberger et alj 2004al) : 



/(logP) 



-(1/2) Z 2 



2vr 



(10) 
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Table 4. List of candidate triple systems. It consists of Sgo systems with P < 400 d, flagged as DMSA/G or A/x systems, with 
no firm evidence so far of being a triple system. In column remar k, 'not resolved' mean s that no visual companion is listed in the 
Fourth Catalog of Interferometric Measurements of Binary Stars dHartkopf et al ] l200lh . 



HIP 


DMSA/G 


P(d) 


e 


Remark 


3362 


7 


2.17 


0.0 


SB2 (dMle + dMle) 


12716 


5 


0.95 





twin SB2: G6 + K0; not resolved 


14273 


7 


2.73 


0.0 


SB 2 


18080 


5 


0.51 


0.07 


W UMa binary 


19248 


7 


2.56 





SB2 


23657 


5 


1.58 


0.06 


SB2 (B5V + B5V); unlikely to be A^i binary because parallax 
is very small 


34630 


7 


20.63 


0.296 


SB2 


45080 


7 


6.74 


0.18 




50097 


5 


3.24 


0.0 




53295 


7 


15.83 


0.31 


not resolved 


70931 


9 


11.82 


0.33 




89601 


7 


4.48 


0.00 


not resolved 


95028 


7 


208.8 


0.37 


not resolved 


98351 


7 


62.88 


0.34 


not resolved 


100142 


7 


2.98 


0.07 


twin SB2: B2V + B2V; not resolved 


104987 


5 


98.82 


0.0 


SB2 & visual binary (G2III + A5V) 


105406 


7 


3.24 


0.02 


SB2 (F8V + K6V); not resolved 


110514 


7 


6.72 


0.178 


twin SB2: F8IV + F7IV 


114639 


7 


3.96 


0.04 


SB2: K1IV + F8V; not resolved 



Weak candidates 



21673 


9 


38.95 


0.15 


not resolved; orbital elements considered 'marginal' 


52444 


7 


51.55 


0.0 


orbital elements based on small number of data points 


68523 


9 


207.36 


0.55 


not resolved; orbital elements uncertain 


77911 


7 


1.26 


0.61 


either triple system as suggested by eccentricity pumping or 
doubtful orbital elements (see text) 




10 20 

parallax(mas) 



Fig. 7. Same as Fig. [5] but for barium stars. 



where z = (logP — 4.8)/s and s — 2.3, periods being ex- 
pressed in days. More precisely, the fraction of binaries with 
periods in the range Ppm (= 1000 d) - -Psb among binaries with 
periods P < Psb may be expressed as 1 — <£>(zpm)/ < £ > (.zsb), 
where 




To obtain the computed fractions 1 — $(zpm) /&(zsb) — 0.36 
and 1 — &(zA[i) /®(zsb) = 0.29 consistent with the observed 
fractions (0.357 ± 0.025 and 0.299 ± 0.024, respectively), it is 
required that zsb = —0.42 (or Psb = 6800 d). The capability 
of the Geneva-Copenhagen survey to detect spectroscopic bi- 
naries up to periods of about 6800 d is not inconsistent with its 
typical radial-velocity sampling (see for instance Fig.|9]l. 

We now evaluate the total number of proper-motion and 
spectroscopic binaries in the Geneva-Copenhagen catalogue. 
In total, 1892 stars out of 9431 with HIP and Tycho-2 entries 

in the Geneva-Copenhagen catalogue have Prob(xy ) < 0.001 
and may thus be considered as bona fide spectroscopic binaries. 
On top of these, 384 stars are DMSA/G or A/x binaries (at the 
0.9999 level) not seen as spectroscopic binaries (Prob(xy ) > 
0.001), for a total of 803 DMSA/G or proper-motion binaries in 
the Geneva-Copenhagen catalogue (or 8.5%). In total, the frac- 
tion of (spectroscopic + astrometric) binaries in the Geneva- 
Copenhagen survey thus amounts to 24.1 ± 0.4% [= (1892 + 
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Table 5. Duplicity of barium stars with parallaxes w > 
5 mas. The flags in the second column correspond to posi- 
tive (or negative) detections with the following methods: (0) 
no radial-velocity variations detected by CORAVEL mon- 
itoring (1) spectroscopic binary (2) astrometric orbit sus- 
pected (3) astrometric orbit constrained from spectroscopic 
orbit (4) A/i binary with Prob(x 2 ) > 0.9999 (5) A/x bi- 
nary with 0.9 < Prob(x 2 ) < 0.9999 (6) accelerated proper 
motion (DM SA/G binary). Flags (0) and (1) are derived 
fromlJorissen et all dl998l).IUdrv et all ( 1998), Nordst r5m et all 
(200 41); flagWfr om Jorissen et al. (2004, 2005); flag (3) 
from Pourbaix & Jorissen (2000), Jancart et al.l ( 120051) ; flags 
(4,5) from the present paper; and flag (6) from the Hipparcos 
Catalogue DMSA. Parallaxes of barium stars flagged as A/i bi- 
naries at a confidence level of 0.9999 (flag 4) have been written 
in bold face. 



HIP 


Flags 


-^orb 




Rem. 






(d) 


(mas) 




944 


26 




7.11 




4422 


_ 


_ 


15.84 




10119 


25 




8.91 


Ba dwarf 


13055 


12346 


2018 


6.54 




20102 


1 


1694 


6.01 




31205 


123 


457 


8.25 




32894 


125 


912 


14.11 




33628 


1256 


> 3700 


6.47 




39426 


05 




6.34 




47267 


1256 




6.73 


SB in Gen-Cop 


49166 






7.90 


50805 


123 


667 


23.42 




53807 







9.54 




56130 






5.46 




56731 


123 


1711 


7.07 




58948 


14 


> 4700 


19.08 




60292 


1235 


3569 


5.78 




62608 


5 




6.20 




65535 


246 




15.73 




69176 






8.75 




71058 


4 




10.63 




76425 


15 


5324 


13.89 




80356 


15 


> 3400 


8.01 




81729 






23.35 




89386 


2 




5.69 




89881 


25 




9.37 




93188 


2 




5.55 




98431 


4 




14.06 




103263 


1 


4606 


6.33 




104732 


1246 


6489 


21.62 




105881 


12 


2378 


8.19 




105969 


13 


878 


16.61 


Ba dwarf 


106306 


125 


2837 


6.31 




107818 






7.53 




112821 


14 


4098 


10.74 




114155 


1 


111 


6.07 




117585 


25 




7.08 





384)/ 9431]. When restricting the sample to stars with paral- 
laxes in excess of 20 mas, this frequency becomes 20.7 ± 0.8% 
[= (364+ 131)/ 2386]. 

The 384 proper-motion binaries not detected as spectro- 
scopic binaries generally have few radial-velocity points (< 5) 
not covering a long time span (At < 2000 d), as can be seen 
on Fig. [9] However, thirty systems do not fall in this category: 
their radial-velocity data had the required properties to con- 
firm their binary nature inferred from the proper-motion data, 
and yet they are not flagged as spectroscopic binaries at the 
level Prob(xy ) < 0.001 in the Geneva-Copenhagen survey. 
These systems are listed in Table |6l under three different cate- 
gories: (i) single stars spuriously flagged as proper-motion bi- 



Table 5. Continued for barium stars with w < 5 mas. Only 
systems for which at least one flag can be set are listed. 



LITTj — 


Mags 


TJ 

-*orb 

(d) 


■uj 
(mas) 


Rem. 


rm — 

558 


1 — r 1 

5 


- 


1 — 7T 

1.65 




21519 


5 


- 


4.55 




21681 


25 


■ 


3.44 




25452 


15 


2652 


3.32 




25547 


25 


- 


3.43 




27767 


5 


- 


0.83 




29740 


123 


1689 


-1.25 




30338 


13 


629 


1.56 




32713 


13 


1768 


0.73 




34143 


15 


7500 


3.65 




34795 


02 


- 


2.44 




36042 


13 


672 


2.36 




38488 


2 


- 


1.29 




42999 


5 


- 


0.99 




47881 


26 


- 


3.56 




50259 


5 


- 


4.38 




51533 


5 


1754 


4.00 




52271 


13 


918 


3.4 




53091 


2 


~ 


2.14 




53717 


15 


1653 


2.30 




61 175 


5 




4.89 




66844 


05 




3.87 




88743 


5 




3.75 




90316 


2 




0.56 




96024 


2 




3.48 




97613 


246 




4.10 




97874 


2 




2.07 




105294 


246 




3.49 


Ba dwarf 


107478 


2 




3.75 




107685 


2 




1.77 




109662 


5 




-0.10 




110108 


15 


1019 


4.59 




113641 


5 




3.17 




117607 


1356 


1294 


4.61 




117829 


5 




4.42 





i i | i i i | i i i 

0.0 <= probaVr < 0.8 




, 0.8 <= probaVr < 0.999 




0.999 <= probaVr 


All 


//T t 




/ j - 






/ 





0.2 0.4 0.6 0.8 

Probfo*) 



Fig. 8. Comparison of the cumulative Prob(x 2 ) (from Eq. [TJ 
for the Geneva-Copenhagen sample of F and G dwarfs, for 
different values of their radial-velocity variability probability 
(probaVr = 1 — Prob(x 2 / ))• The largest fraction of proper- 
motion binaries is indeed found among radial-velocity binaries, 
having probaVr > 0.999. 
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naries because they exhibit perspective accelerations, (ii) con- 
firmed spectroscopic or visual binaries, or (iii) stars suspected 
of being long -period binaries (and seen nearly pole-on, or with 
low-mass companions, or with very eccentric orbits). We dis- 
cuss each category in turn in the remainder of this section. 

The first category corresponds to a purely geometrical ef- 
fect, caused by the perspective diminishing (increasing) the 
proper motion of a close, fastly receding (respectively ap- 
proaching) object. Writing the transverse spatial velocity as 
Vt — fid = \V\ cos (3 and the radial velocity as V r — \V\ sin/3 
(where dis the distance to an object, \V\ — (Vf+Vj?) 1 ' 2 is the 
modulus of the spatial velocity, (3 is the angle between the tan- 
gential and spatial velocities, such that $ = fi) and expressing 
that, for a single star, |V| remains constant over time, yields 



/i = —2 /iV r d~ 



(12) 



This perspective acceleration /i may be detected as a secu- 
lar change of the proper motion, and confused with a proper- 
motion binary. The variation A/i of the proper motion over the 
interval At may be written as: 



|A/t| = \2[iV r d- 1 At (1 - V r d- 1 At) |, 



(13) 



corresponding to the second-order Taylor expansion of the so- 
lution of Eq. [T2] Perspective acceleration should be detectable 
by our x 2 test with a threshold of 18.4206 for stars having 
| A/i| > 4.3 (TA/i, where 0A/j is a typical error of a or 6 com- 
ponent of the proper motion difference between Hipparcos and 
Tycho-2. In practice, the situation is a bit more complicated, 
since 'instantaneous' proper motions (as those appearing in the 
previous formula) cannot be determined. Proper motions are 
always computed from differences of positions extending over 
some time interval, and the effect of perspective acceleration 
on the positions is proportional to /i At 2 . The Tycho-2 proper 
motion is derived by a linear least-square fit on these positions 
spanning At ~ 100 yr: the Tycho-2 proper motion will there- 
fore fall somewhere in the range defined by the 'instantaneous' 
proper motions corresponding to each end of the time interval 
At, so that | /imp — /iTycho-2| < \Afi\, with |A/t| expressed 
byEq.Ha 

There is no need here, however, to go beyond the simple 
approximation expressed by the inequality \2 /t V r d _1 At\ > 
4.3 cta/z (with At ~ 100 yr and 0A/j ~ 1-4 mas y _1 , as suit- 
able for Hipparcos and Tycho-2 measurements), since the pur- 
pose of the whole discussion is just to identify stars showing 
perspective accelerations mistaken for A/i binaries. These stars 
are listed in Table[6](and were already identified as perspective- 
acceleration stars in Ta ble 1.2.3 o f Volume 1 of the Hipparcos 
and Tycho Catalogues; lESAlll997h . For none of the stars listed 
in Table|6]is there a need to go beyond the second-order Taylor 
expansion considered in Eq.[T3l 

All but four stars in Table [6] do not, by far, fulfill the 
criterion for perspective acceleration. An extensive literature 
search for these 26 stars resulted in 8 stars with clear indica- 
tions for being close visual binaries or spectroscopic binaries. 
Interestingly enou gh, two among these 8 systems have low- 
mass (HIP 29860: IVogt et~alll2002h or even brown-dwarf com- 
panions (HIP 29295 = Gliese 229: lNakaiima et alll995l) . They 




N(Vr) 



Fig. 9. Properties of the radial-velocity data for stars detected 
as proper-motion binaries at the 0.9999 level: N(V r ) is the 
number of radial-velocity data points and At their time span 
(in days). Open circles correspond to systems detected as spec- 
troscopic binaries (SB) at the level Prob(x 2 / ) < 0.001, and 
filled triangles to proper-motion binaries with constant radial 
velocity at the level Prob(x 2 / ) > 0.001. It is clear that most 
of the proper-motion binaries not detected as spectroscopic bi- 
naries have either radial-velocity data points that are scarce 
(N(V r ) < 5), or that have a short time span (At < 2000 d). 



constitute further examples of the ability of the A/i method 
to detect binaries with low-mass compan ions, adding to the 
case o f HIP 25647 already discussed by iMakarov & Kaplan 
d2005l) . In the case of HIP 80902 (=HD 150706) which has 
been announced to host a 1 Jupiter-mass exoplanet dUdry et al. 
2003b . one may however doubt that the A/t method is 



able to detect companion s that light (see the discussion of 
Mak arov & Kaplan! 120051 in relation with their Eq. 2). The 
18 stars in Table [6] neither exhibiting perspective accelerations 
nor already known as binaries are flagged candidate binaries 
with possibly low-mass companions, although very low incli- 
nations or very eccentric orbits could be other explanations for 
their non-detection as spectroscopic binaries by the Geneva- 
Copenhagen survey. 

5.4. K and M giants in the solar neighbourhood 

Fama ev et al. (120051) mblished a study of the local kinematics 
of K and M giants, accompanied by a catalogue that contains 
6691 such objects, possessing Hipparcos and Tycho-2 proper 
motions, as well as CORAVEL radial velocities. The sample 
can thus be subject to a similar analysis as done above for the 
Geneva-Copenhagen catalogue. Again, a threshold of 0.999 on 
l-Prob(xv;J to identify stars with variable radial velocities has 
been adopted. The confrontation of spectroscopic and A/i bi- 
naries must, however, be restricted to the 5680 K giants, since 
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Table 6. Stars flagged as DMSA/G or (candidate) A/i-binaries, not detected as spectroscopic binaries in the Geneva-Copenhagen 
survey despite enough radial-velocity data points (N(V r ) > 5) spanning more than 2000 d. For spectroscopic binaries detected 
elsewhere, Psb is the spectroscopic period, if applicable. If A"/i[mas yr _1 ] \V r \ [km s _1 ] ro[mas] > 3.0 10~ 2 mas yr~ 2 , where 
K = 10~ 9 [kpc km -1 s yr _1 ] is a unit conversion factor to get fi in mas yr~ 2 , perspective acceleration is detectable. 



HIP 


Name 


N(V r ) At(V r ) 


Prob( X "(K)) Kfimp \V r \ m 


Psb 


Remark 






(d) 


(mas yr~ 2 ) 


(d) 





Perspective acceleration 



24186 


Kapteyn 


6 


4368 


0.195 


0.54 




54035 


HD 95735 


19 


5487 


0.081 


0.16 




99461 


HD 191408 


9 


2459 


0.939 


0.035 




104217 


61 Cyg B 


32 


6226 


0.723 


0.094 




iinary systems 












12421 


HD 16619 


6 


3940 


0.280 


1 10" 4 




29295 


HD 42581 


15 


5488 


0.003 


5 10~ 4 




29860 


HD 43587 


9 


2953 


0.010 


1 10~ 4 


12309 


60268 


HD 107582 


6 


4676 


0.805 


8 10~ 4 


6864 


67275 


t Boo 


33 


5618 


0.002 


5 10- 4 


3.3 + . 


80902 


HD 150706 


12 


4048 


0.670 


8 10 _s 


264 


102431 


HD 198084 


14 


6625 


0.036 


3 10~ 4 


523 


110618 


v Ind 


16 


3196 


0.612 


1 10~ 3 





visual binary (0.1") 

Gl 229 B brown dwarf 

(Nakaiima et al. 1995b 

Mcomp > 0.34 M Q i lVogt et all 

12055) 

SB JLatham et al1l2002h 
exoplan et in a w id e bi- 
nary ([Pa tience et al. 2002; 
Egaenberger et al. 20 04bT) 
exoplanetary companion 
JUdrv et alJ2003T) 
SB2 <GriffiiiT999l) 
visual binary (0.1"): A3V? + 
F9 



Candidate binaries with low-mass companions (or nearly pole-on systems or very eccentric systems) 



493 


HD 101 


7 


2452 


0.924 


2 10~ 4 






4446 


HD 5357 


6 


4375 


0.128 


1 10~ 5 






8674 


HD 11397 


6 


2206 


0.923 


2 10~ 4 


CH 


subdwarf 



16788 
19335 



HD 22309 
HD 25998 



13 



5247 
2908 



0.045 
0.234 



2 10 

3 10 



44896 


HD 78612 


45 


2567 


0.008 


8 10" 4 


65268 


HD 116316 


7 


2288 


0.571 


8 10~ 5 


66319 


HD 118186 


50 


2438 


0.226 


8 10" 6 


91215 


HD 171706 


6 


2113 


0.685 


5 10 _s 


95703 


HD 183473 


6 


3615 


0.003 


8 10" 5 


98192 


HD 189087 


20 


5194 


0.980 


3 10~ 4 


101726 


HD 196141 


6 


3299 


0.085 


8 10~ 5 


103260 


HD 198829 


9 


3734 


0.015 


1 10 -5 


103859 


HD 200560 


12 


2961 


0.918 


3 10" 4 


108095 


HD 208068 


7 


2644 


0.932 


9 10" 5 


109461 


LHS 3767 


11 


2835 


0.015 


4 10~ 4 


111004 


HD 213122 


6 


2227 


0.111 


2 10~ 7 


112117 


HD 214953 


6 


5429 


0.714 


2 10" 4 



( Sleivvte & Bartkevicius 

[i99Tn ^ 

non SB in lLatham etakl J2002T) 
no low-mass companion 
detected with ELODIE 
(Galland et al. 2005b 



in M giants intrinsic radial-velocity variations (related to enve 
lope pulsations) make the Prob(x 2 / ) cri terion useless to ide n 
tify spectroscopic binaries (see Fig. 5 of iFamaev et al ] l2005l) . 



This section presents an analysis similar to that performed 
for the Geneva-Copenhagen sample of F-G dwarfs (Sect. 15.31 
compare Fig. [10] with Fig. [8j for K giants. The fraction of 
A/i binaries (at the confidence level of 0.9999) among spec- 
troscopic binaries [Prob(xy ) < 0.001] amounts to 5.6% (= 
46/816), raising to 21.8% when adding the DMSA/G sys- 
tems. These fractions both become equal to 53.8 ± 13.8% (= 



7/13) when eliminating those stars with parallaxes smaller than 
20 mas among which proper-motion binaries are not efficiently 
detected. The latter frequency differs by about 1.5er from the 
corresponding values for the Geneva-Copenhagen sample, and 
a possible cause for that difference will be suggested in Sect. [6] 
The large errors on the K-giant frequencies reflect the fact that 
they are derived on a small sample, since not many K giants 
have parallaxes in excess of 20 mas. 

We now evaluate the total number of proper-motion and 
spectroscopic binaries among K giants. In total, 816 stars out 
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0.2 0.4 0.6 0.8 1 10 20 

Prob(x 2 ) parallax(mas) 



Fig. 10. Same a s Fig. [8] but for the K giants from the sample 
of|Famaey et al. (2005). The largest fraction of proper-motion 
binaries is indeed found among radial-velocity binaries, having 
probaVr > 0.999. 



of 5680 have Prob(xy ) < 0.001 and may thus be considered 
as bona fide spectroscopic binaries. On top of these, 82 stars 
are DMSA/G or A/i binaries (at the 0.9999 level) not seen as 
spectroscopic binaries (Prob(xy.) > 0.001), for a total of 178 
DMSA/G or proper-motion binaries. The fraction of (spectro- 
scopic + astrometric) binaries among K giants thus amounts to 
15.8% [= (816 + 82)/ 5680]. This fraction raises to 30.2 ±6.3% 
[= (13 + 3)/ 53 ] when restricting the sample to K giants with 
parallaxes in excess of 20 mas. This frequency again differs by 
1.5cr from the corresponding 20.7 ± 0.8% frequency derived 
for the Geneva-Copenhagen sample. 

Among the 82 proper-motion binaries not detected as spec- 
troscopic binaries, only two (HIP 24450 and HIP 78159) have 
more than 5 radial-velocity measurements spanning more than 
1000 d, thus with the potential to detect a possible long -period 
orbit. Actually, HIP 78159 is a visual binary with 1.8" angu- 
lar separation and a period of 41.56 y (about 15 000 d; see 
Dvorak et alJ [l989. and references therein) which most prob- 
ably corresponds to the proper-motion binary. No independent 
confirmation is available in the literature about the binary na- 
ture of HIP 24450. 



6. Fraction of proper-motion binaries as a 
function of spectral type 

Unbiased estimates of the proportion of binaries among differ- 
ent spectral types are rare, so that the opportunity offered by 
the A/i and DMSA/G binaries is worth investigating. The effi- 
ciency of detecting A/i binaries is affected by the distance to 
a tested object. Among farther objects more will be missed by 
the method due their smaller (angular) proper motions. To ob- 
tain a reliable estimate of the binary fraction, it is necessary to 



Fig. 11. Upper panel: The parallax distribution of stars with 
spectral type G in the Hipparcos catalogue. Lower panel: 
Fraction of stars detected as A/i binaries (at confidence levels 
0.9 - dashed line - and 0.9999 - dotted line) in a given parallax 
bin. Note how the dotted curve reaches its final level of 10% 
for parallaxes larger than 20 mas. At 10 mas, only half of the 
expected systems are detected. 



restrict the sample to the distance (i.e., parallax) region where 
the A/i method attains its full efficiency. Figure [TTl displays 
the detection fraction as a function of parallax for G type stars 
in our complete sample. It can be seen that a plateau in the 
detected binary fraction is reached beyond 20 mas. In the fol- 
lowing analysis of the A/i binary frequency among the various 
spectral types, the sample has therefore been restricted to stars 
with parallaxes larger than 20 mas. 

Table Q lists the fraction of proper-motion binaries among 
the various spectral types. Because A/i and DMSA/G bina- 
ries probe slightly different period ranges (A/i binaries not 
flagged as DMSA/G, as listed in Col. 5, represent the long- 
period tail among proper-motion binaries: see Fig.|U), they have 
been listed separately in Table [7] The importance of restricting 
the sample to stars with parallaxes larger than 20 mas is very 
clearly illustrated by the comparison of Cols. 1 and 2. The frac- 
tion of proper-motion binaries stays basically constant (within 
the error bars) among the various spectral types (Cols. 2 to 5) 
when restricting the sample to parallaxes larger than 20 mas, in 
contrast to the situation prevailing for the full sample (Col. 1). 
Since the proportion of stars with w > 20 mas (Col. 6) de- 
pends upon the spectral type (because stars of different spectral 
types - and hence average luminosities - have different average 
distances in a magnitude-limited sample), the fraction of the 
sample in which proper-motion binaries are efficiently detected 
varies with spectral type, resulting in the variations observed 
in Col. 1 of Table [7] S ince this selecti o n bias was not prop- 
erly recognised by iMakarov & Kaplan! d2005l) . these authors 
concluded that the fraction of long-period binaries is larger 
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Table 7. Fraction of proper-motion binaries of a given type among various spectral classes, normalised by the total number 
of stars in the considered spectral class (labelled ./V tot ). Uncertainties on the listed frequencies are derived from the binomial 
distribution. The large fluctuations observed in Col. (1) are not physical and are caused by the different parallax distributions 
among the different spectral types [see Col. (6) and text]. 



Spectral type 


iVtot 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) Catalogue 


WR-0 


251 


0.012±0.007 


< 0.004 


< 0.004 


< 0.004 


< 0.004 


0.4% 


B 


8664 


0.022±0.001 


0.143±0.132 


0.143 ±0.132 


0.143±0.132 


< 0.0001 


0.8% 


A 


15662 


0.033±0.001 


0.152±0.030 


0.131 ±0.028 


0.083±0.022 


0.069±0.021 


0.9% 


F 


21340 


0.046±0.001 


0.114±0.010 


0.091 ± 0.009 


0.063±0.008 


0.051 ±0.007 


4.3% 


G 


19628 


0.043±0.001 


0.105±0.008 


0.091 ± 0.007 


0.052±0.006 


0.053±0.006 


8.2% 


K 


29349 


0.026±0.001 


0.093±0.007 


0.070 ± 0.006 


0.051±0.005 


0.042±0.005 


6.2% 


M 


4217 


0.029±0.002 


0.098±0.013 


0.077 ±0.012 


0.036±0.008 


0.061±0.010 


12.3% 


All 


103304 


0.0346±0.0005 


0.101±0.004 


0.080 ± 0.004 


0.053±0.003 


0.048±0.003 


5.3% 



Am 


1356 


0.075±0.007 


0.167±0.048 


0.133 ± 0.044 


0.117±0.041 


0.050±0.028 


4.4% 


(a) 


Gen-Cop (F-G V) 


9431 


0.085±0.003 


0.109±0.006 


0.089 ± 0.006 


0.057±0.005 


0.052±0.004 


25.3% 


(b) 


Gliese (G-K V) 


236 


0.106±0.020 


0.109±0.021 


0.087 ±0.018 


0.043±0.013 


0.065±0.016 


97.5% 


(c) 


Kill 


5680 


0.031±0.002 


0.189±0.054 


0.189 ±0.054 


0.132±0.046 


0.057±0.031 


0.9% 


(d) 


Ba 


156 


0.090±0.023 


0.27±0.07* 


0.19 ±0.06* 


0.16±0.06* 


0.11±0.05* 


23.7%* 


(e) 








0.50±0.14** 


0.50 ±0.14** 


0.17±0.11" 


0.33±0.14** 


7.7%** 


(e) 



Column headings: 

(1) A/i (Prob(x 2 ) > 0.9999) + DMSA/G (All Prob), all vj 

(2) A/i (Prob(x 2 ) > 0.9999) + DMSA/G (All Prob), vj > 20 mas 

(3) A/i (Prob(x 2 ) > 0.9999), w > 20 mas 

(4) DMSA/G, All Prob(x 2 ), w > 20 mas 

(5) Only A/i (Prob(x 2 ) > 0.9999), no DMSA/G (All Prob), w > 20 mas 

(6) Fraction of stars in the considered spectral class with w > 20 mas 

* For barium stars, these frequencies refer to a parallax threshold of 5 instead of 20 mas 
** For barium stars, these frequencies refer to a parallax threshold of 10 instead of 20 mas 

Ca talogue reference: 

(a) Renson et al. ( 1991) (b) Nordstrom et al. (2004) (c) Halbwachs et al. (2003) (d) Famaey et al. (2005) (e) Lii et al. ( 1983) 



among solar-type stars than among the other spectral classes 
(as is indeed suggested by Col. 1 of Table |7). Based on the 
above discussion, this conclusion must be discarded, however, 
as shown by the nearly constant frequencies listed in Cols. 2 
to 5 of Table [7] For instance, the fraction of proper- motion bi- 
naries among F-G stars is now about the same for the Geneva- 
Copenhagen sample, for the whole Hipparc os Catalogue and 



forG and K dwarfs from the Gliese catalogue (Hal bwachs et al 
2003). 



If anything, the fraction of proper-motion binaries seems 
(slightly) larger among spectral types associated with high lu- 
minosities, like the upper main sequence (B and A types) and 
the K giants. This effect may again result from a selection bias, 
remembering that proper-motion binaries cannot be detected in 
binaries with components of nearly equal luminosities (since 
in that case, the photocentre will be located close to the centre 
of mass, and there will not be any detectable orbital motion): 
in binaries with a high-luminosity primary star, there is a larger 
luminosity range accessible for the secondary (belonging to the 
lower main sequence), thus increasing the probability of detect- 
ing those systems as proper-motion binaries. 

In the same vein, it should be remarked on Fig. [TT] that 
there is a large fraction of A/i binaries among stars with neg- 
ative parallaxes (at least when the confidence level is relaxed 
from 0.9999 to 0.9), suggesting that the same disturbing ef- 
fect (unrecognised orbital motion?) that pushed the astrometric 
solution into unphysical territory may have led to (slightly) in- 
correct Hipparcos proper motions. Finally, it is worth nothing 
that, for positive parallaxes, the fraction of A/i binaries at the 



0.9 confidence level runs almost parallel to the curve corre- 
sponding to 0.9999, the two curves being offset by 10%, corre- 
sponding to the fraction of false detections expected at the 0.9 
confidence level. This constitutes a further confirmation of the 
validity of the statistical treatment described in Sect. 

Barium stars are an interesting comparison sample since 
they are known to be all binaries (see the discussion in 
Sect. 15.2b . and yet the frequency of A/i binaries listed in 
Col. (1) of Table |7] is only 0.09 ± 0.02 because of the small 
proportion of barium stars with large parallaxes (only 2 stars 
have parallaxes larger than 20 mas). To get meaningful frequen- 
cies, it has thus been necessary to lower the parallax threshold 
adopted for barium stars in Table|7](see also Sect. l5.21 >. Two dif- 
ferent values have been adopted: 5 mas, to ensure consistency 
with the discussion of Sect. 15.21 and 10 mas to get a somewhat 
larger detection efficiency of proper-motion binaries. Although 
significantly larger - as expected - than the frequency of bi- 
naries found in the other spectral classes, this frequency is not 
unity because (i) the 5 and 10 mas parallax thresholds only 
ensure a partial (about 30 to 50%; see Figs.|5]|2l andfTTTi detec- 
tion of the proper-motion binaries; (ii) as already discussed in 
Sect. 15. 31 A/i binaries are restricted to a limited range in orbital 
periods (about 1500-30000 d, as clearly apparent from Fig. |4), 
whereas barium binaries extend beyond that period range. 

Am stars are another interesting class to evaluate the fre- 
quency of A/i binaries, in the framework of a possible rela- 
tionship between this clas s and the barium stars dHakkilal 19891 



Jorissen & Van Eckll2005l) . The arguments in support of such a 



connection are based on the binary nature of stars from both 
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classes, on the similarities of their abundance peculiarities and 
of their kinematical properties. The main argument to reject 
any such evolutionary link, however, is based on the signifi- 
cantly different period ranges observed among Am and barium 
stars. However, the question arises whether long-period bina- 
ries are really lacking among Am stars, or whether this lack 
is a consequence of the difficulty of identifying them from a 
radial-velocity monitoring. It is possible to check this alterna- 
tive by deriving the fraction of A/i binaries amo ng the sample 
of 135 6 Am (and Ap) stars from the catalogue of Renson et al.l 
(fl99l ] laving an HIP entry. That fraction is the same as for A- 
F dwarfs so that there is no hint for more long-period binaries 
among Am stars than among the comparison sample. 

S stars without lin es from the unstable element Tc (also 
called extrinsic S stars; Van Eck et all 19981) are a family of bi- 
nary stars clos ely related to barium stars dJorissen et al.l l 1998: 
Jorissenl 12003). It is therefore of interest to check the frac- 
tion of A/i binaries among them as well. The situation is, 
however, less favourable for extrinsic S stars than for barium 
stars because (i) extrinsic S stars are more luminous (hence 
they are more distant on a verage and thus hav e smaller par- 
allaxes) than barium stars dVan Eck et al.l ll998). and (ii) they 
are mixed with intrinsic S stars (i.e., genuine asymptotic gi- 
ant branch stars), which need not be binaries. Actually, none 
among the extrinsic S stars in the Hipparcos Catalogue has a 
parallax in excess of 10 mas, and it is therefore not surpris- 
ing that no A/i binary is present among them. The only A/i 
binary detected among S stars is HIP 110487 (ir 1 Gru), an in- 
trinsic S star that is k nown to be a visua l binary with a GOV 
companion 2.7" away dProust et all 1981 ). for wh ich an orbital 
period of 6000 yrs is inferred ( Kna pp et al . 1999). With such a 
long orbital period, one may actually wonder whether the vi- 
sual companion is really the one responsible for the A/i binary, 
or whether there could be another, closer companion that plays 
some role as well in the shaping of the complex circumstellar 
environment inferred fro m the su bmm CO lines bvlKnapp et al.l 
1999). IChiuet~ai1 ( l2006l) . and lCruzalebes & Sacutol d20Qi 



The search for A// binaries among long-period variables 
(LPVs) offers an interesting alternative to radial-velocity mon- 
itoring since spectroscopic binaries are very difficult to detect 
among LPVs, due to the intrins ic radial-velocity jitter caused 
by their pulsation djorissenl2003l) . Table[8]reveals a few Afi and 
DMSA/G binaries amo ng the 233 LPVs fr om the Hipparcos 
Catalogue collected by Knapp et"aD d2003h . These candidate 
proper-motion binaries are unusual in the sense that they are 
detected among stars with very low parallaxes. In fact, among 
LPVs with parallaxes in excess of 10 mas (only 3 stars), there 
are no candidate binaries. This casts doubts on the reality of 
the binary detection for stars with parallaxes as small as those 
listed in Table [8] Except for o Cet, there is unfortunately no 
independent confirmation of their duplicity. Parallaxes of very 
red stars such as LPVs are subject to specific instrumental ef- 
fects like chromaticity, causing a displacement of their stellar 
image on the detector (as compared to the image location of a 
blue r star occupying th e sa me position on the sky ), as discussed 
by lKnappet"atl d2003l) and IPlatais etaD d2003l) . It is possible 
that this chromatic effect, if not properly corrected on the older 
plates used by the Tycho-2 consortium to derive the proper mo- 



tion, could spoil its determination and produce spurious A/i 
binaries. 

Finally, to conclude this section, it is interesting to list the 
few candidate A/i binaries found among Wolf-Rayet (WR) 
and O-type sta rs, adding to the sp ectr oscopic and vi s ual bi - 
naries listed by lMason et al.l dl998 ) and lNiemela et al.l d 1999b . 
Individual properties of A/i binaries among WR and O stars, 
down to a confidence level of Prob(x 2 ) > 0.9 9, are listed in 



Table [9] T his list may be com pared to that of lHartkopf et al 



(1999) and lMason et alj d!998l) . who searched for visual bina- 
ries among O and WR stars using speckle interferometry. These 
speckle observations could not confirm the duplicity of any of 
the stars of Table|9] Here as for LPVs, the small parallaxes cast 
doubts on the reality of the candidate binaries. 

7. Conclusions 

The possibility of detecting binaries from the inconsistency be- 
tween short-span (like Hipparcos) and long-sp an (like Tycho- 
2) proper motions (an idea put forward b y IWielenl Il997 : 
Wielenetal1ll998l; iMakarov & Kaplan! 120051) has been inves- 
tigated using test samples of known (spectroscopic) binaries. 
Using Hipparcos, Tycho-2, and S^o data, we show that almost 
all binaries in the period range 1000 - 30000 d (the lower bound 
corresponding to the span of the Hipparcos mission) are in- 
deed detected as 'proper-motion binaries'. In the lower range 
1000 to 2000 d, binaries are best detected from their accel- 
erated proper motions (flagged as DMSA/G in the Hipparcos 
Double and Multiple Systems Annex). Proper-motion binaries 
with 5*89 periods below 1000 d are good candidate triple stars, 
the proper-motion binary involving a component with a longer 
orbital period than the Sga period. A detailed literature search 
indeed confirms the triple nature of 16 systems, and a list of 19 
candidate triple systems is provided. 

The comparison of short-span and long-span proper mo- 
tions thus offers a promising way to detect binaries in ex- 
tended stellar samples, and has a bright future given the var- 
ious upcoming space astrometry missions like SIM and Gaia. 
Proper-motion binaries also fill the gap between spectroscopic 
and speckle/interferometric binaries. When estimating the fre- 
quency of binaries in a given sample, their contribution is not 
negligible: in the Geneva-Copenhagen sample of F and G stars, 
proper-motion binaries not detected as spectroscopic binaries 
add 5.5% to the fraction of spectroscopic binaries, 15.2%, for 
a total of 20.7 ± 0.8%. In the Geneva-Copenhagen sample, 
there are several proper-motion binaries not detected as spec- 
troscopic binaries, despite adequate radial-velocity coverage. 
These binaries may either be pole-on systems, very eccentric 
systems or systems with low-mass companions. 

The frequency of proper-motion binaries in fact appears 
to be the same among all spectral classes (about 10.1 ± 
0.4% in the whole F-M ra nge), contrasting with the claim of 
Makarov & Kaplan! d2005l) that it is higher among solar-type 
stars (because these authors did not account for the selection 
biases). There is instead a hint that slightly more proper- motion 
binaries are observed for systems with luminous primaries (AV 
or Kill), as expected if the light ratio with their companion 
is generally larger than for systems hosting primaries on the 
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Table 8. Candidate A/i binaries at the 0.9999 confidence le vel or DMSA /G binaries among long-period variable stars. The Col. 
labelled tn(rev) lists the parallax recomputed as explained in Knapp et al. (2003). See text for a word of caution about the reality 
of these candidate binaries, inspired by their very small parallaxes. 



HIP 


GCVS 


Prob(^) 


x 2 


DMSA 


ccr(HIP) 


(rev) 












(mas) 


(mas) 


10826 


o Ceti 


1.0000000 


48.66 


5 


7.79 


9.34±1.07 


14042 


THor 


0.9999998 


30.94 


V 


-0.15 


0.69±1.22 


28874 


S Lep 


0.9999836 


22.03 


5 


3.63 




34859 


VX Gem 


1.0000000 


46.72 


5 


1.05 


1.02±1.99 


78872 


ZSco 


0.1404664 


0.30 


G 


-2.47 


1.71± 2.56 


84948 


RS Her 


0.9309712 


5.34 


G 


1.10 


0.80±1.41 


85617 


TW Oph 


1.0000000 


37.73 


V 


3.57 


3.66±1.17 



Table 9. Candidate A/i binaries at the 0.99 confidence level or DMSA/G binaries among Wolf-Rayet and O-type stars. The Col. 
labelled m lists the HIPPARCOS parallax. SB in Col. 'remark' stands for 'spectroscopic binary'. See text for a word of caution 
about the reality of these candidate binaries, inspired by their very small parallaxes. 



HIP 


HD 


Prob( X ^) 


X 2 DMSA Sp. type 


cc(HIP) Remark 










(mas) 



48617 


86161 


0.9921 


9.69 


5 


WR 16 


-1.01 ±0.74 


SB a 10.73 d 


68995 


123008 


0.9963 


11.21 


5 


O 


1.95 ± 1.30 




89769 


168206 


0.9989 


13.62 


5 


WR 113 


1.38 ± 1.32 


SB 6 29.1 d 


99768 


192639 


0.9926 


9.83 


5 


O 


1.22 ±0.64 




114685 


219286 


0.9994 


14.98 


5 


07 


2.04 ± 0.90 





(a) Moffat & Niemela ( 1982) (b) Massev & Niemela ( 1981); Niemela et al. ( 1999) 



lower main sequence. Hence, the photocentric motion is ex- 
pected to be more frequently detectable - all other parameters 
being equal - for the former than for the latter. 

As expected, barium stars exhibit a significantly larger frac- 
tion of proper-motion binaries (50±15% for stars with za > 
10 mas, consistent with 100% when duly accounting for the 
~ 50 % detection efficiency of proper-motion binaries when 
vd ~ 10 mas). There is no hint at a similarly large fraction of 
long-period binaries among Am stars; the suspicion of a link 
between these two classes must thus be definitely abandoned. 
No convincing candidate proper-motion binaries emerged ei- 
ther among WR and O stars, or among long-period variable 
stars. 
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